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We calculate the effective charge for multimagnon infrared 
(IR) absorption assisted by phonons in a perovskite like an- 
tiferromagnet and we compute the spectra for two magnon 
absorption using interacting spin-wave theory. The full set 
of equations for the interacting two magnon problem is pre- 
sented in the random phase approximation for arbitrary total 
momentum of the magnon pair. The spin wave theory results 
fit very well the primary peak of recent measured bands in 
the parent insulating compounds of cuprate superconductors. 
The line shape is explained as being due to the absorption of 
one phonon plus a new quasiparticle excitation of the Heisen- 
berg Hamiltonian that consists off a long lived virtual bound 
state of two magnons (bimagnon). The bimagnon states have 
well defined energy and momentum in a substantial portion 
of the Brillouin zone. The higher energy bands are explained 
as one phonon plus higher multimagnon absorption processes. 
Other possible experiments for observing bimagnons are pro- 
posed. In addition we predict the line shape for the spin one 
system La2Ni04. 
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I. INTRODUCTION 

In the late 50 's Newman and Chrenko presented a pi- 
oneering infrared absorption study in NiO Among 
other results they showed a band at 0.24eV that cor- 
related with the disappearance of antiferromagnetism 
above the Neel temperature and hence of likely magnetic 
origin. This was very puzzling because in the NaCl struc- 
ture, with an O at an inversion center, a direct magnetic 
absorption is not allowed In fact in a typical two 
magnon excitation, like for example, two spin flips on 
adjacent metals atoms, all relaxation of charge in the 
excited state is symmetric with no net dipole moment. 
An explanation of this puzzle was given by Mizuno and 
Koide 1^ who observed that the joint absorption process 
of a phonon and two magnons is allowed since then the 
symmetry of the lattice, after the phonon excitation, is 
eff'ectively lower. To the best of our knowledge no de- 
tailed theory existed of this effect. Note that the Mizuno 
and Koide ||] paper is from the early days of Anderson's 
superexchange theory. 

Also very puzzling data was recently presented by 
Perkins et al |^,^. They measured the absorption in 
many different parent insulating compounds of high-Tc 



superconductors. The data shows a narrow primary peak 
in the charge transfer gap and a set of side bands. Ini- 
tially the narrow peak was associated with an uniden- 
tified exciton. However no exciton is expected in this 
energy range 

The propose of this paper is to present a detailed the- 
ory of phonon-assisted multimagnon infrared absorption. 
We apply the theory to layered insulators and show that 
it explains Perkins et al. data (Sec. fv\) - We estimate the 
coupling constant of light with multimagnon excitations 
(Sec. H) and we calculate the line shape for one phonon 
plus two magnon absorption using interacting spin wave 
theory. For this we need to solve the two-magnon prob- 
lem for arbitrary total momentum. We present the 
full set of equations in the random-phase-approximation 
(RPA) and solve them approximately (Sec. [II). The nar- 
row primary peak reproduced in Fig. ^ is explained in 
terms of a new quasiparticle excitation. It consists of a 
long lived virtual bound state of magnons referred to here 
as a bimagnon. The new state has well defined energy 
and momentum in a substantial portion of the Brillouin 
zone. We also present a prediction for the line shape of 
the spin one system La2Ni04. 

So far complementary information on antiferromag- 
netism in perovskite materials has come from probes like 
neutron scattering and Raman light scattering Our 
results show that IR absorption can be used to study 
magnetic properties of this systems. This is interesting 
because the technique is relatively simple. 

In the conclusions (Sec. 0) we discuss and propose 
other experiments which should be able to detect the 
new excitations. 

A short account of this results was presented elsewhere 



II. THEORY OF PHONON-ASSISTED IR 
ABSORPTION OF MAGNONS 

A. Model Hamiltonian 

We concentrate on the case of a spin 5=1/2 and 2 
dimensional (2-d) material like the Cu-0 layers of the 
cuprates but we indicate the corresponding generaliza- 
tions for a system with larger dimension and/or spin 
like La2Ni04 or NiO. We consider a three-band-Peierls- 



1 



Hubbard model ||T^.|TT| in the presence of an electric field 




Here, cj^ creates a hole with spin a at site i in the Cu 
d^2_y2 or the O Px,y orbital. For simplicity Cu atoms 
are kept fixed and O atoms are allowed to move with 
displacements Uk- These restrictions will be efi^ectively 
removed when we compare with real experiments since 
measured phonon properties will be used. ICki is a spring 
constant tensor and Pi is the momentum canonically 
conjugate to ui, Pph = Ze'^^Uk is the phonon dipo- 
lar moment and Z is the ionic charge of O {Z = —2). 
For electron-lattice coupling, we assume that the nearest- 
neighbor Cu-0 hopping is modified by the 0-ion displace- 
ment Uk as tij = t ± aduk, where 5 is a unit vector in 
the direction of the corresponding Cu-Cu bond, 6 = x, 
y (and z for a 3-d system) and the + (— ) applies if the 
Cu-0 bond shrinks (stretches) with positive 5uk (0-0 
hopping is neglected here), t is given by l/2V3pd a m 
terms of Slater-Koster integrals The site energy 

contains the coupling to the electric field and in addition 
the Cu site energy is assumed to be modulated linearly 
by the displacements of the O ions, 

e,^Ed + pJ2^±)Suk + eEr, (Cu), (2) 

A; 

e,^Ep + eEr, (O). (3) 

In Eq. (|^) the sum extends over the four surrounding 
O ions. The sign takes the value + (— ) if the bond 
becomes longer (shorter) with positive Suk and r; is 
the position of atom i (including displacement). The 
other electronic matrix elements are: Cu-sitc (Ud) and 
0-site {Up) repulsions for Ui, and the nearest-neighbor 
Cu-0 repulsion (Upd)- We define A = Ep — Ed + Upd 
and e = 2{Ep — Ed) + Up. In the following we adopt 
the notation that Cu sites are labeled with i and O 
sites or the corresponding Cu-Cu bond with i + S/2. 
In this notation the position of an O ion is given by 
f'i+s/2 — '2(* + ^/2) +Mi+6/2j with a the lattice constant. 

B. Effective Charges and Dipole Moment Operator 

To calculate the coupling constants of light with one- 
and two-phonon-multimagnon processes we first obtain 
a low energy Hamiltonian as a perturbation expansion 
valid when t << A, e, Ud and when the phonon field and 



the electric field vary slowly with respect to typical gap 
frequencies, 

H = J2 JiE: {ui+s/2})Bi+s/2 + - -BPph. (4) 

i,S 

Here i3i+6/2 — SiSi^s with Si spin operators, Hp^ is 
the phonon Hamiltonian containing spring constants and 
masses for the O ions (M) and Pph- The first term 
in Eq. (^) contains the spin-dependent fourth order cor- 
rection in t whereas fourth, second and zero order spin- 
independent processes are collected in the last two terms. 
To compute J we can use the three center system Cul- 
0-CuR of Fig. I 

See Ref. |] 



FIG. 1. Schematic representation of the cluster used in the 
calculations. Full dots represent Cu's and open dots O's. 
Thick arrows represent the spin, thin short arrows represent 
lattice displacements and thin long arrows represent the di- 
rection of the electric field. We have represented uq in con- 
figuration A. In general its direction is equal to the direction 
of the electric field. 

In order to generalize our result for arbitrary spin S 
we define the saturated ferromagnetic state 

|F) = I'^L = S, Si, = S; Tor = S, S'r, = S) 

and the Neel state, 

|N) — |tol — S,Si, — S; Tor — —S, S'r, = S). 

Here tor^l is the z component and S'l,r is the magnitude 
of the total spin of the ion. The superexchange can be 
calculated as, 

J = -^{(FIVRVRVRVIF) - (mVRVRVRVm) (5) 

where R — (1 — P)/ {Eq — Hq), P projects on the manifold 
with Sh = Sk = S, Ho contains all site-diagonal terms 
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of the Hamiltonian Eq. (||) and V the d — p hybridization 
terms. For details see Ref. In the case of a 3-d 

S=l system hke NiO we can orient the z axis in the 
direction of the Ni-O-Ni bond, then t should be taken 
as the hybridization between (12^2 orbitals and the Pz 
orbital and is given by pda in terms of Slater and Koster 
integrals Also appropriate values of Ud should be 

used The same is valid for La2Ni04 if the crystal 

field splitting between the d^^^^ri and the dx^-y^ orbitals 
is neglected. 

We only need to consider the three configurations 
(A,B,C) of the L-R bond and the electric field. Next 
we Taylor expand J to first order in E and second order 
in {ui+s/2}, 

J = Jo + r]{uL - wr) - E[qiUo + AqA(2uo - - ur)] 
- E[^iuo{ui, - Mr) + ACa(2uo - Ml - '"r)('"l - ur)] 
+ ... (6) 

Here A = 1 for configuration A and A = for configu- 
rations B and C. In each configuration the displacement 
of the central O and the electric field are parallel, i.e. 
E = Ee,Uo = uqE. The direction of e is the same 
as the arrows at the bottom of Fig. |l|. ul and ur are 
only relevant in configuration A. ul ~ uli + ul2 — ul3, 
Mr = — MRi-f-UR2-f ur3. The numbering and the direction 
of the displacements are shown in Fig. Other terms 
quadratic in u's have been neglected, they renormalize 
the spring constants in iJph and give the magnon-two- 
phonon interaction. The first term in Eq. is the su- 
perexchange in absence of the electric and phonon fields. 



Jn = 



I 



(7) 



The remaining quantities are a magnon-phonon coupling 
constant, 
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-(7r + 



(8) 



effective charges associated with one phonon and multi- 
magnon processes, 
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(9) 



2i'' I I 2 2 

f"^ 9 12 11? 

-e-i^/3a,,[^(7r + -) + 7r(x + 77-)']' ^ 



^A2^C/<j t' C//A Ud' 



and effective charges associated with fwo-phonon and 
multimagnon processes. 

t^(3 .Z ,1 2, 8,1 1,2, 

(11) 

2t^l3'^apd, 2 3 4 2 3 



apd is the Cu-0 distance, a/2. Within a point charge esti- 
mation the parameter Papd ~ 2Upd- The dipole moment 
is obtained from Eq. (||) as P = — and using Eq. (||) 
in the relevant configurations. We get up to fourth order 
in t, 

P = Plph + f 2ph + Plph+mag + P2ph+mag (12) 

The first two terms describe conventional one and two 
phonon absorption processes. 

We define 6Bi+s/2 = -Bi+6/2- (-Bi+6/2) and its Fourier 
transform. 



(13) 



The lattice spacing a and h are set to f . The supraindex 
labels the direction of the vector d and N is the number 
of unit cells. In the same way the Fourier transform of 
Ui+s/2 is given by 



u" 



i+S/2- 



(14) 



After Fourier transforming, the dipole moment for one 
phonon and multimagnon processes for an in-plane field 
in the x direction is, 



Ipli+mag 



p ^xp 



(15) 



pS 



and A = 1. For an electric field perpendicular to the 
plane we have, 



Iph+mag 
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E 

pS 



(16) 



In a 3-d system like NiO, A = 1 and we get the analogous 
of Eq. ( [isl ) for the three directions and the sum over S is 
also for the three directions. 

The term proportional to qi is isotropic being present 
in any configuration. Looking at the cluster in Fig. |l| it 
can be understood as a spin dependent correction to the 
charge on Oq. Its physical origin is that fourth order cor- 
rections to the charges involve spin dependent processes. 
For example if the spins in Cul and Cur are parallel 
they cannot both transfer to Oq whereas if they are an- 
tiparallel they can. Fig. ||(a) illustrates a typical process 
efficient in configuration B. 

An alternative way of deriving qi is with the aid of the 
Hellmann-Feynman theorem. We can add to the elec- 
tronic Hamiltonian of the Cul-O-Cur a term eno where 
n-O — '^OT + "-Oi and noa is the occupation number op- 
erator for the central O ion and spin a. Now we can 
calculate the expectation values of the resulting Hamil- 
tonian H{e) in the ground state |e). We have that the 



3 



total charge on the O is given by the following derivative 
of the energy, 



d{e\H{e)\e) , 
|?io|0) = e — le^o- 



(17) 



Since the term proportional to Jo{e) in (e|iJ(e)|e) is the 
magnetic part of the energy its derivative is related to 
the magnetic dependent part of the O charge, 



dJoje) 
' de 



\6=0, 



and since e just renormalizes A we have formally 

dJo 



qi = e 



dA' 



(18) 



(19) 



This expression is useful to obtain qi if one can get Jq 
from a different method. 

See Ref. il 



potential in Cur decreases, and in Cul it increases, cre- 
ating a displacement of charge from left to right that 
contributes to the dipole moment. Fig. ||(b) illustrates 
a typical process. But again this effect is spin depen- 
dent since if the two spins are parallel they cannot both 
transfer to Cur, and hence one gets a spin dependent cor- 
rection to the dipole charge. In a highly covalent material 
like the cuprates one expect the charge phonon effects to 
be quite strong and dominate the effective charges (see 
Sec.|rVA|). 



C. Optical absorption 

The real part of the optical conductivity due to the 
processes described in the previous sections is given by 
the dipole-moment-dipole-moment correlation function. 
In Zubarev's notation [[1(1 



2ttuj 
NVcu 



Im(((Piph+mag; ^Iph+mag))) 



(20) 



Here Vcu is the volume associated with a Cu ion, i.e. 
NVcu is the total volume of the system per Cu-0 layer. 
Now we assume for simplicity that only the ufp's with the 
same 6 and e mix. To zero order in the magnon-phonon 
interaction we can decouple the magnetic system from 
the phonon system. In this approximation the eigen- 
states of the system are products of phonon states times 
magnetic states and by writing the Lehmann representa- 
tion of the Green function one can factor out all phonon 
matrix elements. Using that 



(Oph I ufpuf_p I Oph) 



1 
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FIG. 2. Typical processes contributing to the isotropic (a) 
and anisotropic (b) eflective charges. The meaning of the 
symbols is the same as in Fig. 2. 

The term proportional to in Eq. ( p^ is anisotropic 
being present for an in-plane field only or when there are 
Cu-Cu bonds oriented in the direction of the field. It orig- 
inates from a "charged phonon" like effect |l^ . Consider 
the configuration in which the electric field and the dis- 
placement of Oq are both parallel to the Cul-Cur bond 
(A in Fig. and a phonon in which the O's around 
Cul breathe in and the O's around Cur breathe out. 
(We don't need to consider zero momentum phonons to 
couple to light since it is the total momentum, magnons 
plus phonons which has to add to zero.) The Madelung 



with |Oph) the phonon vacuum and wfp the phonon fre- 
quency. We get. 



MVcu 



16A2(7isin2(2|^)sin2(Si) + {4.XqASm^i^) - qif 



OJ 



\p 



(21) 



Here t^Hp is the frequency of the u'^p and u^p phonons and 
ll}±p is the frequency of the u^p and Uyp phonons. t^Hp 
can be associated with the frequency of Cu-0 stretch- 
ing mode phonons and Lu±p with that of Cu-0 bending 
mode phonons. The supraindex in the Green functions 
indicates that the poles should be shifted by that amount. 
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III. THE TWO MAGNON PROBLEM IN 
INTERACTING SPIN WAVE THEORY 

To compute the magnon-magnon Green functions we 
use interacting spin-wave theory with a Holstein- 
Primakoff transformation. On the A sublattice we put 



'2S'(1- 



2S 



)bi 



and on the B sublattice, 



2S 



(22a) 



St= b] 



'25(1-|^) 



b]bi 



-S + b\bi 



(22b) 



where bi is a Boson operator. 

Now the Hamiltonian can be formally expanded in 
powers of 1/S. From now on we adopt a more conven- 
tional notation and drop the subindex in Jg- We define 
rii = b\bi. The Hamiltonian reads H = i^Ncci + Ho + Hi. 
With, 



Ho 



rij 



(23) 
(24) 



<^J> 
J 



Hi = -jJ2 ^i^j I'^W'^i + + '^•C'] (25) 



<ij> 



<^3> 



< ij > indicates that nearest neighbor pairs are counted 
once in the sum, z is the coordination number, £^Neei 
is the classical Neel energy, Hq is the linear spin wave 
Hamiltonian and Hi is the spin-wave-spin-wave interac- 
tion. Notice that the Hamiltonian is invariant under the 
exchange of the sublattices and so we don't need to dis- 
tinguish between them. Accordingly we work in the non- 
magnetic Brillouin zone. 

The non- interacting part, Ho is diagonalized by the 
Bogoliubov transformation. 



A = "fc^fc ^ Vkb-k 



Qk = Ukbk - Vkb^_^. 
where bk is the Fourier transform of bi and 



(26) 



Uk 



Vk 



1 + ^fc 

2iLJk 



-sig(7fe 



(27) 



2uok 



and 



Ik = - y^cos(fea), 



(28) 
(29) 



Now we can normal order the Hamiltonian with respect 
to the non-interacting spin-wave ground state. This is 
equivalent to writing the Hamiltonian as a Hartree-Fock 
part plus fluctuations which we latter treat in the RPA. 
After normal ordering the Hamiltonian can be written 
as, 

H = -EnccI + Esw + E HkQkQk + Kcs, (30) 

k 

where Egw = jJSNzdl + C/2S), Ek = E^,uok. E,,, = 
zSJ[l + C/25'), C is the Oguchi correction, 



k 



(31) 



For z = 4 as in the layered materials, C, ~ 0.158. Vros = 
v}is + V^res contains the normal ordered product of the 
interacting part i.e. 



-J E] • niJij : 

<ij> 

J 



(32) 



(33) 



<^3> 



The first term originates from the Ising part of the in- 
teraction. In the Ising limit it is easy to check that the 
effect of this term is to shift the energy of two nearest- 
neighbors spin flips from the noninteracting value 4 J to 
3J. In the Heisenberg limit and in the case of the Ra- 
man line shape as an example, the non interacting line 
shape has a peak at 4J and this term shifts it to close 
to 3J. This term is the more important one to get the 
correct line shape. The second term is a correction to 
the exchange due to the kinematic interaction. 

The next step is to put V-^cs in terms of the spin-wave 
operators Eq. (p6|). Then we evaluate the normal order 
and we do the RPA approximation. This consists in this 
case, in keeping only those terms in Vros which create or 
destroy a pair of magnons. With this we get. 



2N 



2 - 3 - 4)ri234Qlg^Q3Q4 (34) 



1234 
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1,2,... stands for fci,fc2,- 
the Appendix 
We define 



and the vertex are given in 



p is the total momentum of a magnon pair and q^, 
the relative momentum. In the following we make indis- 
criminate use of the property that gpqiq2 does not change 
when <7i —^ —gi- 
lt is convenient for latter use to define the following 
auxiliary functions 



a. 



the following "different sublattice" form factors, 
/pq = a+gCOs{qs), 

fpq = ttpq C0S((75), 

f =/?+ 

J pq r^pq ' 

and "same sublattice" form factors. 



(36a) 
(36b) 



(37) 



hpq = /3+qCOs((75), 

<q=/3pqSin(g5), (38) 

= a+ 

PQ pq 5 

pq pq' 

The /'s have the property that /pq+Tr = ^fpqj with 
TT = (tt, TT, ...). If we Fourier transform them in q we 
see that they are different from zero in different sublat- 
tices. For the /I's, hpq^^^ = hpq and they are different 
from zero on the same sublattice. They will allow us 
to classify the Green functions in "different sublattices" 
Green functions and "same sublattice" Green function 
where the relevant coordinate is the distance between the 
spin operators. For example in the case of the operator 
Eq. (^3|) since S joints different sublattices we need the 
former. 

By replacing Vies by l^^J'^ in the Hamiltonian Eq. ( |30| ) 
the RPA equation of motion for the Green functions can 
be computed by standard methods ||l^. No further ap- 
proximations are needed to obtain the Eqs. (^),(45). We 

get 



1 



^t^" ~ ]y 2^ gpqlT-pqn-pq 

q 

^lJ.v = "77 Qpqq'hpqhpqi 



(40) 
(41) 



qq' 



where /i, = 1^5 for I = 1,2 and — I for Z = 3,4. 
Analogously the different sublattice Green functions are 
given by, 



'^t^t'— ^ A^ypqJpqJpq 

q 



^'"^ Y 9pqq' fpqfpq' 



(42) 
(43) 



qq' 



It is easy to see that similar zero order Green functions 
with an / and an h factor in the kernel vanish. Since 



the vertex [Eq. (A3)] do not mix same sublattice with 



different sublattice Green functions the equations for the 
interacting Green functions separate in two blocks. From 
Eq. (|3^) we get the following for same sublattice Green 
functions, 



2 

-sin{—)K'l^l]K2Sv 



(44) 



and for different sublattices, 



G^. = G^S-^lEtC^ + cos(^)G(°3)]Gi,. 



s 



+ [^lpG% + \Y.^os[^)G';:>,]G,^}. (45) 

5 

Note that the Green functions depend on the total mo- 
mentum p which appears also as a parameter in the equa- 
tions. 

We can expand also our two magnon Green function 
in Eq. ( pT] ) in powers of l/S" to lowest order we get. 



(0) 



9pq.q. - '^"'"+/"-" g°q. + ^ E ^P.. rpq, qffpqq,, ((^^^-p ; ^^^p)) = ^ 



(39) 



and the vertex 



where o° = — - Ei„_„ , 

is given in the Appendix lA^. Since the vertex is a sum 
of separable potentials the equations can be solved in 
integral form. We define same sublattice Green functions. 



[Gi,i, + cos(^)(G3i. -f Gi,3) + cos2(^)G33] (46) 

This neglects small Oguchi type corrections to the pref- 
actor ]T7[ | which slightly renormalizes the intensity scales 
in the results that follow. First we will solve the equa- 
tions in some limiting case in which one can handle the 
equations analytically. 
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A. 2-d Raman 

It is useful to see how some well known results are 
recovered in this formalism. In Appendix |^ we solve the 
equations for the case = Py The case of p = is of 
particular interest because it corresponds to the Raman 
line shape. According to Ref. for the case of a layered 
antiferromagnet the scattered intensity is proportional to 
the imaginary part of the following Green function, 



g2 

— Gr 

TT 



N{{B^ 



^11, p= 



(47) 



Gfj, is defined in the Appendix M and, in this case 



where, 



g: 



,<i(0) 
22,p=0 



and G?(°)^o = L^'^ 



N ^ 



f. 



N ^ {Ujqf UJ - 2Eg 



(48) 



and fq = cosqx — cosqy. With this equivalences Eq. (B^) 
becomes equivalent to the well known expression p^|, 



Gr = 



i(2) 



4[L(0)L(2) 



(L(i)) 



1 + i[m + L(2)) + J±[mm - (m) 



(49) 



Note that the fact that we use a Holstein-Primakoff 
representation rather than a Dyson-Maleev one (as in 
Ref. ||l^) did not affect this result. This is because the 
effect of the exchange part of the interaction which is 
different in the two formalisms cancels out in this case 
due to symmetry. In Fig. |^ we show with a solid line the 
imaginary part of Gr which gives the Raman line shape. 
The position of the maximum is at 



Ef^"^^^ = 1.46£;,n = 3.38J 



(50) 



as is well know ||l7|. Since the line shape lies at high 
energies (w > E'm) we can get a good approximation for 
[Eq. (p3)] by applying the HE approximation of 



^11 



Appendix g. We have G+'+' ~ 0^2^ ~ G 



where in this case, g: 
we get. 



Gp 



d{0) 

++,p=o 



22 

i(L("). 



G 



d(0) 

++,p=o 



f + JG' 



d(0) 



•-^11 > 
.2L(i)+L(2)) and 



(51) 



,p=0 



Here the distinction mention in Appendix ^ between the 
HE approximation done in the vertex alone or in the 
vertex and in the operator does not apply because the 
operator has been explicitly constrained to have the cor- 
rect symmetry. The corresponding line shape is given by 
the short dashed line in Fig. ^ for the Raman case. The 
position of the maximum is at IA6E^ as before and only 
the intensity decrease a bit. A popular and more drastic 



approximation consists in putting Uk 
this case one gets for the Raman case. 



Gr = 



i(o) 



1 + jm ' 



1, Vk = 0. In 



(52) 



This line shape is shown with the long dashed line. The 
peak shifts appreciably and is at 1.48-Bni- This approx- 
imation is much worse and we will not use it in the fol- 
lowing. 




FIG. 3. Raman line shape as given by ImGa according to 
Eq. (^9|) (solid line), in the high energy (HE) approximation 
(short dashed line) and taking Ufe = 1, i;*, = (long dashed 
line). 

It is interesting to compare the present approach with 
the pioneering work of Elliott and Thorpe . They 

assume a Neel state as starting point and neglect the 
Oguchi correction. This gives in this case jlSj a similar 
line shape but peaking at 2.71J. The present theory 
assumes a spin-wave ground state and takes the Oguchi 
correction into account. Although, at least qualitatively 
a similar result is obtain for p = in the two theories, 
this changes dramatically for p 7^ 0. If one tries to extend 
the Elliott and Thorpe theory for that case, nonphysical 
poles appear in the Green functions at energies Eip_^q^ — 
Ei„ „ indicating that a better ground state must be 

2P— qi o o 

used as done here. 



B. IR case 

To calculate the IR line shape we need to evaluate the 
Green function given by Eq. (^|). The narrow exper- 
imental line shape for the cuprates suggests the occur- 
rence of a sharp resonance for some values of p. This can 
not occur close to p = (0, 0) or (tt, tt) because there the 
continua of two magnon excitations extends from zero 
to 2ii^„i and the imaginary part of the non interacting 



7 



Green function is relatively large at the energies of typ- 
ical two-magnon excitations. This changes at different 
values of p, in particular for p — (tt, 0) there is a gap in 
the spectrum from to and the imaginary part of the 
non interacting Green function is very small where the 
resonance occurs. This case is analyzed in detail next. 



1. Bimagnon at p = (tt, 0) 

The relevant Green function for this problem is 
[Eq. m] 



{{SB'- 



Nn 



(53) 



In the Appendix ^ we solve Eq. ( ^ ) for this momentum 
in any dimension. Gixix is given by Eq. (C2). In Fig. ^ 
we show with a solid line the imaginary part of this Green 
function. A sharp resonance occurs indicating that a vir- 
tual bound state (bimagnon) is formed. The maximum 
is at 



E^J^^^ = 1.179Sn 



2.731J 



(54) 



^'[nG^2x2x This repre- 



We show also ImG'J'!^^^., ImG^'l^ja;' 
sents the density of states of the continua of two magnon 
excitations in which the bimagnon can decay. As mention 
before it starts at and is very small at the position 
of the bimagnon pole. This explains the long lifetime of 



the bimagnon. We see that 



d(Q) 



xlx 



G 



d{0) 
2x2x 



^d(O) 
'-'lx2x- 




0.4 0.8 1.2 



FIG. 4. ImGi^io; at p = (vr, 0) in the full theory of Eq. (|Cg) 
(solid line), and almost indistinguishable, in the HE approx- 
imation (doted long dashed line). We also show (very close 
to each other) IniG^^^^ (short dashed line), ImGj^ja; (long 
dashed line), ImGj^ja; (doted short dashed hne). 



This is because at high energies Vk is quite small [see 
Eqs. (36a),(^),(|4^)] and all of them are very close to 
G+l+x defined in Appendix. ^. We can then apply the 
HE approximation of Appendix. ^ For simplicity here 
we use the HE approximation both in the vertex and 
in the operator since in this case there is no low energy 
spurious spectral weight i.e. we neglect G', 



-\-X — X 7 



G'_ 



in Eq. ( p6a ). We get the familiar RPA like form. 



r' 



'~^+x+x 



1 + 2JG 



(0) ■ 



(55) 



The imaginary part of this is shown with doted long 
dashed line in Fig. ^. Again we see that the HE approxi- 
mation is very accurate. In particular the peak does not 
shift and there is only a small decrease in the intensity. 



2. Bimagnon for general p 

To obtain the total line shape we have to integrate 
the contributions from the whole Brillouin zone so we 
need the Green function at all values of p. By the same 
arguments as before we can use the HE approximation 
of Appendix We expect this to work well because of 
the following reasons, i) We expect that the total line 
shape will be determined by the sharp excitations that 
occur close to momentum (tt, 0) and we have seen that in 
this region the approximation does extremely good, ii) 
Far from this region the approximation is also good as 
we have already shown for the Raman case. 

To analyze the following results it is useful to look first 
at the non-interacting {S = oo) case (Fig. H). The regions 
where the imaginary part is small are good candidates for 
having narrow resonances in the interacting case. 



20 r 



(TT.OjO 
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20 r 




FIG. 5. lm{{SB^p-SBp)) in the noninteracting (5 = oo) 
case as a function of p, from (0,0) to (■7r,0) (a); from (7r,0) to 
(tt, tt) (b); and from (0,0) to (7r,7r) (c). The regions in which 
the Une has not been draw have zero imaginary part. 

In Fig. ^ we show the imaginary part of the Green 
function from Eq. (^6|) as a fiinction of the total momen- 
tum p in the high energy approximation in the vertex but 
not in the operator (see Appendix ^). The bimagnon is 
only well defined close to (7r,0). It disperses upwards on 
going towards (0, 0) [Fig. ^(a)] and downwards on going 
towards (tt, tt) [Fig, ^(b)]. This is shown more clearly in 
the inset of Fig. |o| This indicates that (tt, 0) is a saddle 
point and hence it should give a Van Hove singularity 
when integrated over p [Eq. (|2l|)]. Because of that the 
position of the peak in the final integrated line shape is 
the same as for the (tt, 0) bimagnon. 

The peak at lower energy is a real bound state. Its 
position, the intensity, and even the existence are not 
very reliable because the position is beyond the range 
of applicability of the HE approximation. If the HE ap- 
proximation is done in the vertex and in the operator this 
peak gets a much larger spectral weight and contributes 
spurious intensity to the line shape (see next Section). 



-En,Im{{,SBi,:,Si3;» 



(1.0)0 




2 



i5„,Im((^-Bi,:,)_B;» 
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FIG. 6. lm{{SBlj,- 5Bp)) for S = 1/2 in the HE approxima- 
tion in the vertex as a function of p, from (0,0) to (7r,0) (a); 
from (7r,0) to (7r,7r) (b); and from (0,0) to (tt, tt) (c). In the 
case of real bound states a small imaginary part (0.02 Em) 
has been added in the denominator of the Green function. 

In Fig. ^we show the imaginary part of the Green func- 
tion as a function of the total momentum p for the case of 
S = 1. This will correspond to a material like La2Ni04. 
The main bimagnon peak gets overdamped in this case 
because the pole shifts up to the region of larger imagi- 
nary part for the noninteracting Green function (Fig. a). 
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field polarized perpendicular to the plane only phonons 
perpendicular to the Cu-0 bond contribute. A quick es- 
timate indicates that the absorption in this direction is 
roughly a factor of 8 smaller than the in-plane contribu- 
tion. From the experimental side the anisotropy seems 
to be larger. One should be aware that the cuprates are 
in a regime where covalency is not small with respect 
to typical gap energies and hence a perturbation in t is 
helpful to identify the important processes and discuss 
trends, but quantitative estimations are to be taken with 
care. The same problem has been pointed out for the 
computation of J Q . As higher orders in t are included 
we expect that the anisotropic contributions grow with 
respect to the isotropic ones. For example the charged 
phonon effects of Fig. ||(b) can become very efficient if the 
second hole forms a Zhang- Rice singlet with the hole al- 
ready present in Cul since that process involves a much 
smaller gap. Note also that the larger the order in t 
the longer the range of the processes that contribute to 
the anisotropic charges whereas only local processes con- 
tribute to the isotropic charge. Longer range processes 
have in general very large form factors. These effects 
should give a stronger anisotropy in accordance with the 
experiments. This suggest that the process involving uj\\ 
are dominant. We enforce that in the following by taking 
qi = 0. 

B. IR line shape 




FIG. 7. lm{{5B'Lj,;SB^)) for 5 = 1 in the HE approxima- 
tion in the vertex as a function of p, from (0,0) to (7r,0) (a); 
from (7r,0) to {n,n) (b); and from (0,0) to {n,n) (c). In the 
case of real bound states a small imaginary part (0.02 _Bni) 
has been added in the denominator of the Green function. 



IV. COMPARISON WITH EXPERIMENTS 

A. Anisotropy 

We first discuss the anisotropy because this helps to 
fix the values of the effective charges. For an electric 



The line shape for phonon-assisted multimagnon ab- 
sorption is given by Eq. ( pl|) . Since the dispersion of 
the phonons is very small |2l| ] with respect to the line 
width we can neglect it and take Einstein phonons. It is 
convenient to define the following function 

I{uj > 0) = -16 X 

^sin2(^)(sin2(^) +sin2(|))Im((<5i?!Lp;<5i?p) (56) 
p 

and I{uj < 0) = 0. Although the line shape is different 
for each cuprate I{uj) is the same. The line shape is given 

by 



where 



(57) 



(58) 



A/VcuW|| 

The absorption coefficient is obtained assuming weak ab- 
sorption as 

47r 

(59) 

ei 

with ei the real part of the dielectric constant We 
also define 
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ao = — ^ctq. (60) 

In Fig. ^ we show the function /(w) for S" = 1/2 in the 
HE approximation, done in the vertex alone and in both 
the vertex and the operator. We see that at HE the two 
give very similar results but at lower energy the latter 
have a much larger feature coming from the real bound 
states. Most of this spectral weight is spurious and can 
be easily eliminated just by neglecting the contributions 
from the real bound states as was done in Ref. For 
comparison we show I{uj) for other values of S. 




FIG. 8. for different values of S in the HE approx- 

imation done in the vertex alone (full line). In the case of 
S — 1/2 we show also the result for the HE approximation 
done in the vertex and in the operator (dashed line) . A small 
imaginary part (0.016-Em) has been added in the denominator 
of the Green function. 

To choose notice that because of the Van Hove sin- 
gularity the position of the maximum in the line shape 
is determined by the position of the bimagnon peak at 
(^,0) [Eq. H)] and is given by c^m^^o) + E^^^^oy Be- 
cause of that, to optimize the position of the maximum, 
it is convenient to fix the Einstein phonon frequency at 

t^ll = W||(^,0)- 

The dashed curve in Fig. ^ shows the theoretical and 
experimental line shape in the HE approximation (vertex 
only) for Sr2Cu02Cl2. Here q) has not been mea- 
sured so we took it as 5% smaller than the measured IR 
frequency ||2^ of the E„ stretching mode as suggesting by 
comparing with La2Cu04. Notice that these differences 
are insignificant anyway. The fitting for the primary peak 
is quite good. The improved theory (HE approximation 
in the vertex alone) has a slightly narrower line shape 
which makes the fit not as good as for the more approx- 
imate line shape of Ref. Still the fit is quite good 
and this is surprising because such a good fit, especially 
for the width was not possible within RPA in the Raman 



case ||2^. In fact the experimental Raman line shape is 
much broader than the theoretical prediction shown in 
Fig. ||. This suggest that RPA is not so accurate in the 
case of S* = 1/2. The fact that we obtain a reasonable 
fit in our case can be partially reconciled with the rela- 
tively bad performance of RPA in the Raman case by the 
fact that a structure that is artificially broadened around 
p = (tt, 0) by a factor of 2 (but still much narrower than 
the integrated line shape) does not change significantly 
the final result, i.e. the final width in the line shape is 
not so sensitive to the width of the bimagnon close to 
p — (tt, 0) so this experiment does not prove or disprove 
the accuracy of RPA for a particular momentum. 

From the position of the maximum we found J — 
0.107ey. An alternative way to estimate J is from the 
position of the maximum in the Raman peak, Eq. (|5^). 
From the measurements of Ref. [|4j we found J = 
O.lOSeV. In this estimation one should again be aware 
that the accuracy of RPA for a 5 = 1/2 system has been 
questioned |2|]. 

We keep for comparison in Fig's. P,|l^ the contribution 
from the real bound states although their position, inten- 
sity and even their existence are not very reliable because 
they are beyond the applicability of the HE approxima- 
tion. 




Photon Energy [eV] 



FIG. 9. Experimental data from Ref. (solid line) and 
theoretical line shape for two-magnon absorption (dashed 
line) in Sr2Cu02Cl2. The dashed doted line is the contri- 
bution to the line shape from the bimagnon at p = (7r,0). 
Parameters are cj|| = 0.061eV J = 0.107eV, ao = 3.72cm"^ 

The data shows also tiny structures above the primary 
peck. These structures can also be due to adsorbates p^] 
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or, as the same authors suggest, they can be phonon side 
bands. In fact the two-phonon plus bimagnon contribu- 
tion is expected at -E^^^g^ plus two phonon frequencies. 
Its spectral weight is much smaller than the one-phonon 
plus bimagnon. For example the isotropic charge is of 
order ^qi and hence the absorption is of order 

("i + 6/2) II p . p U 

^ W-' Iph+magi -I Iph+mag//) 



2ph+magi -^2ph+mag 



where we have taken a factor ^j"* ~ 1. This is roughly 
three orders of magnitude smaller than the previous con- 
tributions and can explain the tiny structures reported 
for Sr2Cu02Cl2. The shift from the primary peak is ~ 
0.04 eV and 0.07eV in good agreement with typical val- 
ues of LUj_ and . Note that only the presence of a sharp 
bimagnon state would make those processes observable. 

In Fig. |o|we show the fit for La2Cu04. Here we found 
Jo = 0.121eV. The corresponding value from the position 
of the maximum in the observed Raman line shape [ ^J23| 
is J = O.llSeV. The same comment as before applies for 
the accuracy of this RPA estimates. The superexchange 
was estimated also from a study of the moments of the 
Raman line shape with the result J ~ 0.128eV. 

The phonon is very anomalous in orthorombic 
La2Cu04 since it splits due to anharmonicities [Q, its 
partner being at =0.06eV at room temperature. Pre- 
sumably this produces the shoulder observed at lower 
energies in the experiments (Fig. |l^), although the dis- 
tance to the primary is larger than expected which may 
be because other phonons are involved. This feature 
was assigned to direct two-magnon absorption [QH] made 
weakly allowed by the lower lattice symmetry according 
to the results of Ref. [^. However we found that the 
dipole moment for this process is directed in the direc- 
tion bisecting an angle made by the Cu-O-Cu bond ||2^] 
and hence can only contribute for a field perpendicular to 
the plane. It is possible that a soft mode of the distorted 
structure gives rise to this effect since this will gener- 
ate an enhanced intensity at lower frequencies according 
to Eq. (^l]). although more theoretical and experimen- 
tal work is needed to clarify this point. Note that the 
shift between the primary peak and the shoulder is in 
the range of phonon energies (~ 0.04 eV). 



^ 200 
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FIG. 10. Experimental data from Ref. ^ (solid line) and 
theoretical line shape for two-magnon absorption (dashed 
line) in La2Cu04. The dashed doted line is the contribution 
to the line shape from the bimagnon at p = (tt, 0) Parameters 
are = O.OSOeV @, J = 0.121eV, ao = 2.88cm-i 



C. Spectral weights 

The measured spectral weight for the primary peck is 
determined by the parameter ao. With Eq. ( |60| ) we can 
converted that to cro- For Sr2Cu02Cl2 we have e = 5.5 
p2[ and we get ctq ~ 0.023r2~^cm~^ and for La2Cu04 e ~ 
6.0 and ctq = O.OlSrj-^cm"!. We can use Eq. (|5|) 
to estimate qa- For La2Cu04, Vcu = 95 x lO^^^cnr 
and we approximate M with an O mass (a reduced mass 
would be more appropriate). We get qa/g = 0.082. For 
Sr2Cu02Cl2 with Vcu = 123 x IQ-^^cui^ we get qa/c = 
0.088. Examining the expression for the effective charge 
we find that is of order qa/c ~ — ~ 0.1 in very good 
agreement with the experimental values. One should take 
into account however that the previous estimates neglect 
the weight in the side bands so that the observed spectral 
weight is in reality larger than the one estimated with 
^a/g ~ 0.1. This leaves room for the effects discussed in 



Sec. IV A 



D. Prediction for La2Ni04 



In this case structural parameters are similar to the 
previous compound, Vwi = 94 x lO'^'^cm^ and ei = 5.4 
however the ratio J/A is an order of magnitude smaller 
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Ul. So we expect qa 



2JUpd 



0.01 which makes 



the absorption two orders of magnitude smaller than in 
the previous compounds. With this value of qx we get 
ctq = 3.7 X 10~''r2~^cm~^. We have predicted the hue 
shape taking the phonon frequency from neutron scat- 
tering experiments (a;|| = 0.066cV [^) and the superex- 
change form Raman experiments (J = O.OSOeV [|l|). 
The expect absorption is plotted in Fig. |ll|. As explained 
before there are no sharp features in the S = 1 case. 

After this work has been completed we send the pre- 
dicted line shape to J. Perkins and in response we 
received from them the experimental data plotted in 
Fig. The agreement is remarkably good for the line 
shape and our rough estimation of the oscillator strength 
happened to be of the correct order of magnitude. Com- 
parison with the previous case illustrates quite obviously 
the remarkably difference between the S* = 1 and the 
5=1/2 system the latter showing the strong side bands 
which require an explanation beyond RPA |27||. 
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FIG. 11. Theoretical line shape for two-magnon absorp- 
tion in La2Ni04 (dashed line). The dashed doted line is 
the contribution to the line shape from the bimagnon at 
p = (7r,0) Parameters are c^n = 0.066eV Q, J = O.OSOeV, 
ao = 6.0 X 10~^cm~^. The dots are unpublish experimental 
data [Ei. 



V. DISCUSSION AND CONCLUSIONS 

The experiments [|o| show also very strong side bands. 
These are expected as the effect of quantum fluctuations 
corresponding to the creation of an arbitrary even num- 
ber of magnons. This is in principle also true for the 
Raman case and in fact a structure attributed to a four 



magnon process was observed in Raman at a similar 
energy to the first side band (note that phonon energies 
are almost negligible here). 

In the Ising limit one can roughly estimate the posi- 
tion of these side bands. The minimum energy to excite 4 
magnons (plaquette configuration) is CL;4niag = 4 J which 
correspond to the first side band (see Fig's. p|,p^). The 
next side band is attributed to a process in which four 
spins are flipped in a column which has an excitation en- 
ergy of Ci;4mag = 5 J in the Ising limit. In general one can 
show that side bands are expected at integer values of J 
as observed experimentally. We have also computed the 
side bands by exact diagonalization of a small cluster . 
The exact result confirms the Green function calculation 
and show the side bands corresponding to higher multi- 
magnon process. The relative weight of the side bands 
seem to be smaller than in the experiments presumably 
because of finite site effects or the presence of other pro- 
cess in the magnetic Hamiltonian as in Raman |^ , ^ . 
Exact diagonalization calculations of the Raman spectra 
Q show also this kind of side bands. 

As the temperature is raised one expects that the struc- 
ture broadens in a way similar to the Raman line shape 
and at the same time a hot phonon band should appear at 
^{tt^o) ~'^\\ corresponding to a process in which a phonon 
of the bath is absorbed and a bimagnon is emitted. It 
would be interesting if this effect can be seen experimen- 
tally. 

It is interesting to point out that one can in princi- 
ple sec the bimagnons for a particular momentum with- 
out the assistance of the phonon in neutron scattering 
experiments p5[ ] or high resolution electron energy loss 
(HREELS) experiments. In the latter for energy losses 
small with respect to the energy of the incident electrons 
and in specular scattering geometry, dipole selection rules 
are valid and one should also see the analogue of the 
phonon- assisted IR line shape. Preliminary HREELS re- 
sults in Sr2Cu02Cl2 show a structure which agrees 
very well with Perkins et al. data. Other possibility to 
see this excitations without the assistance of the phonon 
would be to add impurities. In this way the zero total 
momentum selection rule will not apply any more and 
the narrow (tt, 0) bimagnon should give some signature 
at iJ^'^Qj . The impurity should also break inversion sym- 
metry (which to same extent is always does) in order 
to make the two magnon excitations IR active and they 
should not introduce doping otherwise the doping depen- 
dent mid-lR band would mask everything. Possible 
candidates would be Nd2-a; Pr^, CUO4 or similar combi- 
nations with other rare earths. Since the Nd is on top of 
a Cu the O is not any more in a center of inversion. 

Another interesting experiment would be to repeat 
Perkins et al. experiment in an insulator with in-plane 
anisotropy. Since an electric field in the x direction cou- 
ples with bimagnons at (tt, 0) and one in the y direction 
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with bimagnons at (0,7r), if this two directions are not 
equivalent the primary peak should split in a twined crys- 
tal. 

We do not believe that this mid-IR excitations are re- 
lated to the doping dependent mid-IR bands observed in 
the cuprates for which a completely different mechanism 
has been proposed p8| , ^. For example for Nd2Cu04_j^ 
in the notation of Ref. [ [40[ the band labeled J at 0.16eV 
has been quite undoubtedly shown to be due to polaron 
formation 0,|ll|,|l|,|2| by Calvani et al. [||. The struc- 
ture labeled / at 0.76eV is the analogue of the mid-IR 
band reported by Uchida et al. for La2-xSra;Cu04 at 
0.5eV and can be explained by a purely electronic mech- 
anism |3^j39t] . Only the much weaker structure labeled 
K at 0.34eV may be related to this excitations. 

In this work we have computed effective coupling con- 
stants of light with multimagnon excitations assisted by 
phonons and the line shape of the primary peak. Our 
results explain recent measured absorption bands in the 
mid IR of parent cuprate superconductors and show that 
this technique proves high energy magnetic excitations. 
This is very interesting because IR spectroscopy is a tech- 
nique intrinsically more accurate and with a much better 
signal to noise ratio than other techniques available like 
Raman or neutron scattering. We have demonstrated the 
existence of very sharp virtual bound states of magnons 
in spin 1/2 systems at momentum (tt, 0) and show the 
shape of the excitations at different momentum. To the 
best of our knowledge is the first time that the RPA 
equations for the two magnon problem are presented for 
arbitrary total momentum. 
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In term of total momentum and relative momentum vari- 
ables the vertex in Eq. (^) can be rewritten as 



P'31'32 ^ / . \-'pqiJpq2 ' pqi pq2) ' nipJpqiJi 



^ 1 

pq2i 



(5i = 1.2 



pqiq2 



(A3) 



2 /\JpqiJpq2 ' 2 

d 

+ \s^v^^(!^)hl^q^Kq2 + 91 ^ 92] 
where we used the definitions in Eqs. (p 



pq2' 



APPENDIX B: SOLUTION OF THE TWO 
MAGNON PROBLEM FOR 2D AND = py 

In this case the Green functions for fixed p are in- 
variant with respect to the exchange <-> Qy. As is 
done in the theory of Raman scattering we can define 
Gfi, = 2{Gix,i'x — Gix,i'y) and a similar definition for 
the noninteracting Green functions. In this case Eq. ( [iF 
reduce to a 2 x 2 problem, 



W 



- E 

2 ^ 

l" = l,2 



^11" ^l"V 



(Bl) 
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APPENDIX A: VERTEX FUNCTIONS 

The symmetrized vertex in Eq. (^) are given by F 

rii + r-L with 

_ II Jz , , , 

r 1234 = - — n2+l{UlV2V:iU4 + ViU2U:iV4) 



which can be solved for 



Gfi 



(B2) 



(Jll + 2\y'll '-'22 ^ ['-'12 ) 



+"-^22 ) + ^22 



APPENDIX C: SOLUTION OF THE TWO 
MAGNON PROBLEM FOR p = (tt.O.O. ) 

In this case the Green function for the x direction 
decouples from the rest in Eq. (^). This is because 
= 0, G'^]i,g = with S ^ X. To see this notice 
that in the definition Eq. (^2|), the kernel changes sign 
under the replacement q ^ q + (tt, 0, 0..). So we get, 



I" xV : 



(CI) 



' = 1,2 



we can solve for 
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Gix1x{tt,Q,0...) 

^(0) , ^(0) 
'-'Ixlx "I" '-'2x2: 



(C2) 



(0) ^(0) 



APPENDIX D: HIGH ENERGY 
APPROXIMATION 

When the hne shape is positioned at sufficiently high 
energies (w > -Em) it is a good app roxim ation to neglect 
the f's respect to the u's in Eqs. ( 36a -38). In principle 
this can be done everywhere in the expressions for the 
Green functions. However it is convenient to do this in 
two steps. First we define 



to the definitions Eqs. (|3q), (|37D,(43), all other contribu- 
tions are negligible. This means that not only the vertex 
or equivalently the interaction part of the Hamiltonian 
(Kos^^) is approximated by (T^o^^'*) but also the opera- 
tor is approximated by ^B'^ . We call this the HE 
approximation in the vertex and in the operator. Al- 
though this is perfectly consistent at high energies this 
approach produces a large spurious contribution at low 
energies. To see this consider the case of p = 0. We can 
define 



SB^ = SB^ - SBl 



(D7) 
(D8) 



±<5 



rl5 _|_ r25 



(Dl) 



and we do the HE approximation in the vertex, i.e. we 
neglect all contributions that involve ■y's in Eq. ([A^). We 
denote with primes the quantities computed in the HE 
approximation, i.e. 



-2J V f+'' f+'' 

S 



(D2) 



r ^ =0 

PQ1Q2 



Now we rederive Eq. ( p^ ) in the HE approximation. 



(D3) 



We can treat (5 as a matrix or vector index of dimension 
d and put this as 



(0) w 

+ '^+1'- 



(D4) 



Notice that G'n, G'22, G'^^ are matrixes G31 are 

vectors and G33 is a scalar. Solving the dxd problem for 
G'_^i we get 



G'^i = (1 + 2 JG 



(0) \-Iq{0) 



+1 ■ 



(D5) 



and replacing in Eq. (D4) we get also G'^,. We also have 
the exact transformation 



G'n — '^iG'+ 



G'_ 



2G' 



G'i3 — 7;{G'+3 + G'_j). 



(D6a) 
(D6b) 
(D6c) 



When calculating the Green function for the two magnon 
operator, Eq. one would be tempted to make the 
same approximation and keep only G^_|_ since according 



Now we find, 



((<5i?-<5i?iJ)) = i ( {{SBf,;SB^)) 
+ {{SB',;dB^,)) 



{{SB^;SB^))). (D9) 



Only the first term contributes to the line shape, i.e 
the p = contribution to the IR line shape is identical 
to the Raman line shape. This follows from the fact 
that Bq is proportional to the Hamiltonian and hence it 
commutes with it. To enforce that in spin- wave theory, 
we would like Bq to commute with H order by order in 
1/5. However if we do the high energy approximation in 
the operator of Eq. (^6|), we miss some terms of order S 
in Bq and then we find that to order S^, [SBp,Ho] ^ 0. 
This gives spurious scattering in the s channel. We can 
avoid that by doing the high energy approximation in the 
vertex alone. In this case G[^, G'^^i and G33 are kept in 
(P) and computed with Eqs. (p4l),(p5|) and (p|). 



Eq, 

It is instructive to check analytically that in this way the 
non-interacting {{SBq;SBq))'^ — 0. 
In the 2-d case we find, 



(1 + 2JG 



with 



(0) N-l 



1 + 2 jg: 
-2jg: 



(0) 
(0) 



-2JG 
f 2JG 



(0) 



x+y+ 
(0) 

X + X+ / 

(DIO) 



A++ = l + 2J(Gi°|,+ + G(V) 

+ 4j2[Git+G(V-(GiV)']- (Dll) 
so we get for G'_^^_^.j. 

I 9 7r^(o) ^(0) _ ^^(0) 

'~'+x+x ' ^-^ r-^+x+x'-'+y+y \'~'+x+y) J 



r' 



(D12) 



This reduce to Eq. ( p5[ ) for p — (tt, 0) using the properties 
of Appendix 
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